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Abstract:  We  study  the  onset  of  synchronization  in  a  network  of  N 

delay-coupled  stochastic  fiber  ring  lasers  with  respect  to  various  parameters 
when  the  coupling  power  is  weak.  In  particular,  for  groups  of  three  or  more 
ring  lasers  mutually  coupled  to  a  central  hub  laser,  we  demonstrate  a  robust 
tendency  toward  out-of-phase  (achronal)  synchronization  between  the  N  —  1 
outer  lasers  and  the  single  inner  laser.  In  contrast  to  the  achronal  synchro¬ 
nization,  we  find  the  outer  lasers  synchronize  with  zero-lag  (isochronal) 
with  respect  to  each  other,  thus  forming  a  set  of  N  —  1  coherent  fiber  lasers. 
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1.  Introduction 

The  dynamics  of  delay-coupled  systems  is  a  rich  and  varied  area  of  study  covering  a  wide 
range  of  topics  from  biological  systems  [1],  to  chemical  reactions  [2],  to  physical  systems  such 
as  laser  systems  [3,4]  and  electrical  circuits  [5].  Understanding  the  behavior  of  such  dynamical 
systems  generally  involves  solving  systems  of  delay  differential  equations,  and  these  systems 
often  exhibit  chaotic  behavior  and  other  complicated  phenomenon  such  as  noise-induced  oscil¬ 
lations.  For  example,  delay-coupled  nonlinear  oscillators  exhibit  complex  phenomena  such  as 
generalized,  phase,  and  lag  synchronization  [6,7]  in  the  chaotic  regimes. 

Systems  of  delay-coupled  oscillators  often  exhibit  a  time  lag  between  individual  nodes,  with 
a  leading  time  series  followed  by  a  lagging  one.  That  is,  lagged  phase  synchronization  oc¬ 
curs  when  the  time  series  are  correlated  but  locked  in  phase  at  a  value  other  than  zero.  Such 
lagged  synchronization  is  called  achronal  synchronization,  and  has  been  studied  for  many 
delay-coupled  dynamical  oscillators,  such  as  mutually  coupled  chaotic  circuits  [8],  and  Ku- 
ramoto  and  Stuart-Landau  oscillators  [9].  In  particular,  achronal  synchronization  is  of  special 
interest  to  those  studying  coupled  arrays  of  nonlinear  optical  systems,  such  as  semiconductor 
or  fiber  ring  lasers  [10].  This  is  due  to  the  fact  that  in  general  coupled  lasing  systems,  the  syn¬ 
chronized  state  with  zero-lag  is  unstable.  Due  to  the  injection  of  noise,  the  unstable  achronal 
state  may  exhibit  spontaneous  switching  between  leading  and  lagging  nodes  [11]. 

There  has  been  substantial  progress  in  understanding  delay-coupled  laser  systems  in  recent 
years.  For  systems  of  delay-coupled  semiconductor  lasers,  achronal  synchronization  has  been 
demonstrated  experimentally  [11]  and  studied  theoretically  [12].  For  erbium  doped  fiber  ring 
lasers  (EDFRLs)  chaotic  behavior  has  been  observed  and  modeled,  demonstrating  a  rich  variety 
of  dynamical  states  [13-15].  Chaotic  synchronization  for  coupled  EDFRLs  has  also  been  ob¬ 
served  and  studied  in  depth  [16, 17],  and  delay-coupled  fiber  ring  lasers  have  also  demonstrated 
properties  of  achronal  synchronization  and  noise-induced  synchronization  [18]. 

Previous  work  has  also  considered  noise-induced  synchronization  in  EDFRLs  connected  via 
passive  coupling  delay  lines.  For  two  coupled  lasers  of  this  type,  both  modeling  and  experi¬ 
mentation  have  demonstrated  a  complex  “leader-follower”  phenomenon  of  dynamic  achronal 
synchronization  between  the  two  lasers,  where  the  offset  of  the  phase  is  precisely  equal  to  the 
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delay  time  introduced  by  the  passive  coupling  line  [19-21].  In  the  case  of  the  two  coupled  fiber 
ring  lasers,  the  achronal  synchronization  is  observed  to  be  non-stationary  in  both  theory  and 
experiment,  since  leading  and  lagging  lasers  switch  roles  over  the  course  of  long  time  observa¬ 
tions.  Such  behavior  is  an  obstruction  to  long  time  synchronization  in  phase  due  to  its  random 
switching  in  an  achronal  state. 

Our  current  work  generalizes  the  EDFRL  model  for  ring  lasers  to  consider  larger  arrays 
than  have  been  modeled  previously.  Specifically,  we  generate  networks  consisting  of  three  or 
more  ring  lasers  attached  to  one  another  through  a  central  “hub”  laser.  In  [22],  it  was  shown  for 
semiconductor  laser  networks  modeled  using  the  Lang-Kobayashi  laser  formulation  [4],  chaotic 
synchronization  occurs  in  the  star  coupling  formation  if  the  number  of  lasers,  N,  is  sufficiently 
large.  Motivated  by  this  result,  we  examined  similar  coupling  arrangements  for  EDFRL  lasers. 
In  contrast  to  the  chaotic  semiconductor  result,  we  show  that  there  is  very  weak  dependence 
on  N.  However,  we  observe  a  uniform  scaling  of  synchronization  onset  as  a  function  of  weak 
coupling.  The  scaling  appears  to  be  independent  of  the  number  of  lasers,  N. 


2.  Mathematical  model 

Using  the  fiber  ring  model  originally  derived  in  [14],  with  modifications  made  in  [23,24]  and 
[19],  we  consider  an  array  of  N  ring  lasers  coupled  via  /V  —  1  injection  lines  comprised  of  a 
passive  single  mode  optical  fiber,  a  splitter,  and  a  variable  attenuator.  The  model  is  characterized 
by  the  total  population  inversion,  of  each  i'h  laser  (averaged  over  the  length  of  the  fiber 

amplifier)  and  the  electric  field  Ej(t)  in  each  laser.  Two  delay  times  occur  in  the  model:  T/<  is 
the  cavity  round-trip  time  in  the  ring,  and  xd  is  the  delay  in  the  coupling  between  lasers.  These 
delay  terms  correspond  to  the  travel  time  of  light  in  various  components  of  the  system,  and 
thus  resolve  the  spatial  dependencies  in  the  system  normally  modeled  by  the  Maxwell-Bloch 
equations.  The  equations  for  the  delay-coupled  dynamics,  including  noise  terms,  are  as  follows: 


where 


Ej(t) 

dWi 

dt 


Rex p  [r(l  -  iaj)Wj(t)  +  iAfl Efdb{t)  +  §,■(*) 


q-  1  -  Wj(t)  -  £fV)  '  {exp  [2  TWj(t)\  -  1 }  , 


1  N 

Ej(t  -  xR)  +  —  £  Bj  kKkEk[t  -  xd). 

k= t 


(1) 

(2) 

(3) 


Ej(t)  denotes  the  complex  envelope  of  the  electric  field  in  laser  j,  measured  at  a  given  reference 
point  inside  of  the  cavity.  Etdb  is  the  feedback  for  laser  j  and  includes  optical  feedback  within 
laser  j  and  optical  coupling  with  other  lasers.  Time  is  dimensionless,  measured  in  units  of  the 
decay  time  of  the  atomic  transition,  yj | .  The  active  medium  is  characterized  by  the  dimensionless 
detuning  ctj  between  the  transition  and  lasing  frequencies,  and  by  the  dimensionless  gain  T  = 
^aL„N().  Here  a  is  the  material  gain,  La  the  length  of  the  active  fiber,  and  No  the  population 
inversion  at  transparency.  Each  ring  cavity  is  characterized  by  its  return  coefficient  R,  which 
represents  the  fraction  of  light  remaining  in  the  cavity  after  one  round-trip,  and  the  average 
phase  change  A <j>  =  2nnLp/X  due  to  propagation  of  light  with  wavelength  A  along  the  passive 
fiber  of  length  Lp  and  index  of  refraction  n.  The  energy  input  for  each  laser  is  given  by  the 
pump  parameter  q.  The  normalization  term  bj  is  the  number  of  lasers  connected  to  laser  j. 

The  coupling  between  lasers  is  through  the  field,  and  described  by  the  connection  matrix,  B. 
Matrix  B  has  components  Bj  k  =  1,  if  node  j  is  directly  connected  to  node  k,  and  zero  otherwise. 
For  the  star  topology,  B  is  symmetric  since  the  coupling  is  assumed  to  be  mutual,  and  given  by 
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the  NxN  matrix: 


(0  1  1  ...  1\ 
1  0  0  ...  0 

1  0  0  ...  0 


\1  0  0  ...  0/ 


(4) 


By  altering  B  in  various  ways,  a  variety  of  different  connection  schemes  for  coupled  laser  sys¬ 
tems  can  be  considered.  The  i %j  terms  model  stimulated  emission  and  are  mean  zero  Gaussian 
noise  applied  to  the  system  with  standard  deviation  D.  Coupling  between  the  lasers  is  charac¬ 
terized  by  the  coupling  strength  Kj,  the  ratio  of  the  power  of  light  in  the  injection  line  to  the 
power  of  light  in  the  source  ring,  and  only  symmetric  coupling  is  considered  at  the  present  time. 


3.  Numerical  methods 

The  model  system  is  simulated  using  an  explicit  time-discretization  and  integration  scheme. 
Following  the  work  of  [19],  the  round-trip  travel  time  Tr  is  calculated  from  the  length  of  the 
active  and  passive  fibers, 

tR  =  {La+Lp)ny\\/c,  (5) 

where  c  is  the  the  speed  of  light.  For  the  purposes  of  discretization,  the  ring  cavity  is  divided 
into  Ns  spatial  elements,  giving  us  a  time  discretization  At  =  Tr/Ns. 

For  each  time-step  of  integration,  the  new  population  inversion  Wj(t  +  At)  is  approximated 
using  a  second-order  modified  Euler  method.  New  values  of  Ej  and  the  feedback  are  computed 
directly  from  the  formulas  1  and  3.  The  raw  data  is  saved  in  a  truncated  time  series  from 
which  the  intensity  Ij  =  \  Ej  | 2  is  computed  and  then  passed  through  a  low -pass  filter.  For  these 
simulations,  a  second  order  Butterworth  filter  from  Matlab  was  used,  and  set  with  a  threshold  of 
125  MHz,  in  keeping  with  experiments  that  use  a  bandwidth  photo-detector  set  at  that  frequency 
[19,24], 

The  results  below  are  generated  using  random  initial  data  for  Ej,  and  initializing  the  inversion 
near  a  steady  state.  Simulations  are  run  for  a  minimum  of  3000-5000  round-trips  in  the  ring 
so  transients  are  removed.  The  parameter  values  used  for  all  simulations  listed  in  Table  1  are 
based  on  previous  experimental  and  modeling  efforts  as  reported  in  [19]. 


4.  Results 


We  explore  the  synchronization  properties  of  arrays  of  fiber  lasers  where  the  structure  of  the 
network  is  in  a  star  configuration.  The  amplitude  cross  correlation  of  any  two  lasers  labeled  a 
and  b  over  a  time  series  of  length  K,  is  defined  by 

Cab  =  —  E  Va{tj  ~  Ts)  —  ( Ia )]  [h(tj)  —  (4)]>  (6) 

Oa<5b  E 


where  Ts  is  a  time  shift  in  the  respective  time  series  of  the  intensities  la  j,  of  two  lasers,  and  <Jaj, 
are  the  standard  deviations  of  the  two  time  series.  We  expect  the  time  series  to  be  offset  due 
to  the  delay  in  communication  based  on  generalized  synchronization  analysis  given  in  [25]  for 
semiconductor  lasers.  Similarly,  we  define  the  phase  coherence  of  two  nodes  as. 


where 


Rab 


K 


j=  i 


EO/ib.j  —  0a  (fy  Ts)  (ph  (tj  ) 


(7) 

(8) 
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Table  1 .  Summary  of  parameters  used  in  the  laser  coupling  model.  The  detuning  param¬ 
eter  of  the  outer  lasers  was  uniform,  aout  =  .0352,  while  a,,,  =  .0202  was  picked  for  the 
inner  laser.  For  simplicity,  the  length  of  the  fibers  (and  thus  the  delays)  were  fixed  for  all 
experiments,  as  were  the  coupling  strengths. 


Parameter 

Value 

Units 

Description 

R 

0.4 

output  coupler  return  coefficient 

a 

2.03  x  10“23 

2 

m 

material  gain  coefficient 

La 

15 

m 

length  of  active  fiber 

Lp 

27 

m 

length  of  passive  fiber 

No 

1020 

m~3 

transparency  inversion 

r 

0.0152 

dimensionless  gain 

O^in 

0.0202 

detuning  factor  (inner  laser) 

0C()iit 

0.0352 

detuning  factor  (all  outer  lasers) 

n 

1.44 

index  of  refraction 

A 

1.55  x  10~6 

m 

wavelength 

A(j) 

1.58  x  10s 

average  phase  change 

D 

0.02 

standard  deviation  of  noise 

q 

40-200 

pump  parameter 

V\\ 

100 

s-1 

population  decay  rate 

zR 

201.6  x  10^9 

s 

cavity  round-trip  time 

L/ 

45  x 10^9 

s 

delay  time  between  lasers 

K 

0-0.014 

coupling  strength 

is  the  phase  difference  between  the  yth  points  of  the  intensities  when  Laser  a  leads  Laser  b.  C 
ranges  between  -1  (negatively  correlated)  and  1  (perfectly  correlated),  while  R  ranges  between 
0  (no  phase  synchronization)  and  1  (complete  phase  synchronization).  Here  the  phase  angle  <j> 
is  defined  as  the  arctangent  of  the  ratio  of  the  real  and  imaginary  parts  of  the  Hilbert  transform 
of  the  intensity  time  series,  following  the  convention  found  in  [19]. 

It  is  inconvenient  to  write  all  the  the  various  Cab  and  Rai,  terms  for  more  than  a  few  lasers, 
and  so  we  introduce  Cout  and  C;„  as  the  mean  value  of  the  cross  correlation  of  all  outer  lasers 
with  one  another,  and  all  outer  lasers  with  the  inner  hub,  respectively.  Analogous  quantities  are 
also  defined  on  R.  To  measure  synchronization  lag  between  between  the  central  laser  and  the 
outer  lasers,  we  examine  which  values  of  time  shift,  zs,  maximize  the  cross  correlation  Cout  and 
On. 

Figure  1  shows  the  mean  cross  correlation  and  phase  coherence  as  a  function  of  zs.  The 
results  were  generated  with  coupling  coefficient  K  =  0.009  and  a  network  of  five  lasers,  four 
on  the  outside  and  one  central  hub.  In  this  simulation  the  four  outer  lasers  are  detuned  from  the 
inner  laser,  with  each  laser  having  the  identical  detuning  factor  a.j  =  .0352,  while  the  detuning 
factor  for  the  inner  laser  is  a,j  =  .0202.  Note  that  when  considering  C,„,  as  in  figure  1,  with 
respect  to  zs  we  adopt  the  convention  that  the  inner  laser  time  series  is  the  one  being  shifted 
backward,  thus  C;„  here  can  be  viewed  as  the  cross  correlation  when  the  inner  laser  is  leading. 
Note  that  in  Figure  1,  the  outer  lasers  have  an  average  peak  cross  and  phase  correlation  at 
zs  =  0,  meaning  that  the  outer  lasers  synchronize  with  each  other.  Concurrently  the  outer  lasers 
synchronize  with  the  inner  hub  at  Zs  =  Zj.  This  is  consistent  with  the  case  of  two  connected 
fiber  lasers,  where  the  synchronization  is  offset  by  one  factor  of  the  delay  time  [11, 19,25]. 

It  is  also  worth  noting  that  an  additional  local  maximum  of  Cout  and  Rout  occurs  at  zs  =  2zj. 
This  agrees  with  intuition,  since  the  total  time  it  takes  for  a  signal  from  one  outer  laser  to 
reach  another  is  exactly  two  delay  cycles.  From  repeated  numerical  experiments  illustrating  the 
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Fig.  1.  The  cross  correlation  and  phase  coherence  for  IV  =  5  fiber  ring  lasers,  as  a  function 
of  a  shift  Ts  in  the  time  series,  for  the  star  configuration.  In  this  experiment,  Xj  =  45ns,  and 
it  is  clear  that  Cout  peaks  at  zero  while  Cm  peaks  at  Tj. 


behavior  described  above,  we  can  conclude  that  the  outer  lasers  tend  to  synchronize  with  one 
another  precisely,  while  the  outer  lasers  tend  to  be  offset  by  a  factor  of  one  delay  with  the  inner 
laser. 

Figure  2  illustrates  the  properties  discussed  above  for  a  simulation  of  5  lasers  in  the  star 
formation.  In  this  simulation,  all  5  lasers  are  given  random  initial  data,  and  the  outer  lasers  are 
identically  detuned  from  the  inner  laser.  The  central  laser,  in  blue,  is  offset  by  Tj,  and  its  peaks 
line  up  with  the  outer  lasers,  all  of  which  are  synchronized. 
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Fig.  2.  A  time  history  of  the  intensity  for  a  numerical  experiment  with  N  =  5  and  K  =  .01. 

Here  the  inner  “hub”  laser  (in  blue)  is  shifted  to  the  right  by  one  delay  because,  in  this 
experiment,  the  outer  lasers  are  being  led  by  the  inner  hub.  Lasers  2-4  are  shifted  up  or 
down  by  a  small  amount  for  illustrative  purposes. 

We  next  explore  how  adding  additional  outer  lasers  affects  the  overall  synchronization  of  the 
system.  Starting  with  N  =  3,  we  run  simulations  with  identical  initial  conditions  and  detuning 
parameters  and  examine  the  final  state  as  a  function  of  K.  Figure  3  shows  how  increasing  the 
coupling  strength  causes  the  lasers  to  go  from  an  unsynchronized  to  a  synchronized  state.  Since 
the  system  (1-3)  behaves  like  a  two  laser  system  when  the  outer  lasers  are  synchronized,  the 
initial  conditions  used  in  generating  figure  3  are  given  by  the  final  state  of  a  two  laser  simulation 
(which  is  simulated  long  enough  for  transient  effects  to  pass).  In  each  experiment,  all  the  outer 
lasers  are  given  the  same  exact  initial  condition  and  prehistory  as  determined  by  one  of  the 
lasers  in  the  two  laser  simulation,  and  the  inner  laser  is  given  the  prehistory  associated  with 
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the  other  laser.  The  error  bars  shown  are  derived  from  doing  batches  of  20  runs  with  different 
random  number  seeds  for  the  noise  terms.  Here,  the  values  of  Cou,  are  taken  with  zero  lag,  and 
were  computed  over  a  time  series  of  ten  round-trips  (Tr).  Note,  as  more  lasers  are  added,  the 
scaling  behavior  of  the  cross  correlation  does  not  change,  but  the  run  to  run  variance  decreases, 
as  evidenced  by  the  smaller  error  bars  on  the  cross-correlation  of  the  outer  lasers.  We  can 
conclude  that  adding  lasers  appears  to  diminish  the  variance  due  to  noise,  but  does  not  enhance 
the  overall  synchronization  scaling  as  a  function  of  K  and  N.  Larger  values  of  N  were  run  (up  to 
N  =  200),  with  no  noticeable  improvement  on  the  overall  synchronization  (for  fixed  k),  though 
these  are  not  shown  in  Figure  3. 


Fig.  3.  Here  we  demonstrate  the  growth  of  the  cross-correlation  for  N  =  3,  N  =  9  and 
N  =  20  lasers  as  a  function  of  K.  Here  Cout  is  recorded  with  Ts  =  0  over  ten  cavity  round 
trips.  Note  that  the  smallest  values  of  K  represent  a  nearly  uncoupled  system. 


Simulations  were  also  run  for  randomly  distributed  detuning  parameters  for  the  outer  lasers 
for  comparison  with  crowd  synchrony  results  obtained  in  [26].  However,  enhanced  synchro¬ 
nization  was  not  observed  for  large  numbers  of  lasers  in  this  regime.  Although  several  factors 
could  be  attributed  to  the  difference  in  results,  it  should  be  noted  that  the  lasers  modeled  here 
exhibit  noise-induced  oscillations,  while  those  in  [26]  exhibit  chaotic  oscillations.  These  obser¬ 
vations  will  be  explored  in-depth  in  a  coming  article. 

5.  Weak  pumping  limit 

As  a  variation,  we  follow  the  work  in  [27]  that  in  the  star  coupling  configuration,  outer  lasers 
tend  to  synchronize,  even  if  the  pump  rate  of  the  inner  laser  is  very  low.  In  fact,  in  certain 
configurations,  the  “weak  link  synchronization”  actually  enhances  synchronization  properties 
when  the  oscillators  are  chaotic.  Although  the  coupling  considered  in  [27]  is  different  than 
what  is  considered  here,  their  coupling  is  evanescent  and  nearest  neighbor,  while  we  focus  on 
injection  line  intensity  coupling,  we  would  like  to  see  if  some  of  these  effects  persist  in  our 
model.  In  Figure  4,  we  turn  the  pump  constant  q  to  be  40  for  the  center  laser,  while  q  =  200 
for  all  outer  lasers,  and  observe  that  the  lasers  still  synchronize  for  small  values  of  K.  In  this 
regime,  the  outer  lasers  synchronize  with  one  another,  even  though  the  amplitudes  of  the  outer 
and  inner  laser(s)  are  not  correlated. 

6.  Conclusion 

We  considered  a  star  coupling  configuration  of  N  coupled  fiber  ring  lasers  in  the  weak  coupling 
limit.  The  model  we  considered  replaced  the  spatial  component  of  the  ring  laser  with  a  time 
delay.  In  addition,  delay  was  included  in  coupling  fibers  used  for  communication.  Using  this 
model,  we  have  shown  synchronization  in  both  amplitude  and  phase  of  the  measured  inten¬ 
sities  can  occur  by  measuring  the  cross  correlation  between  inner  hub  and  outer  lasers.  When 
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Fig.  4.  Achronal  synchronization  for  the  weak  pumping  case.  Here  the  outer  lasers  have 
q  =  200  while  the  inner  have  q  =  40.  Here  K  =  .01.  The  inner  laser  here  is  in  blue  and  is 
blown  up  in  Figure  4b.  Here  Cou,  ~  .901. 


examining  cross  correlation  between  outer  lasers  only,  we  find  that  synchronization  occurs  with 
zero-lag.  This  is  consistent  with  the  results  reported  in  [27],  where  the  outer  lasers  of  evanes- 
cently  coupled  fibers  synchronize  in  the  absence  of  delay.  In  contrast  with  the  outer  lasers,  the 
cross  correlation  between  the  hub  and  outer  lasers  exhibits  lag  synchrony  with  the  lag  time 
equal  to  the  coupling  delay  between  lasers.  That  is,  as  a  group,  the  outer  leaders  lead  or  lag  the 
hub  with  approximately  the  same  waveform  on  average. 

In  the  weak  pumping  limit  as  the  pump  of  hub  approaches  threshold,  we  observe  similar 
results.  As  a  function  of  the  number  of  lasers,  N,  we  do  not  see  a  threshold  for  synchronization 
onset,  however,  the  type  of  coupling  considered  here  is  different  from  that  studied  in  Lang- 
Kobayashi  laser  models  in  [26],  so  perhaps  this  result  is  to  be  expected.  Physically,  we  under¬ 
stand  the  onset  of  synchronization  by  equating  the  array  to  a  mutually  delay  coupled  system 
of  three  lasers,  as  reported  and  analyzed  in  [25].  In  that  work  the  end  lasers  synchronize  while 
they  appear  uncorrelated  with  the  inner,  or  hub,  laser.  The  idea  is  that  the  hub  laser  is  in  general¬ 
ized  synchronization  with  the  outer  lasers,  through  a  nonlinear  function  which  is  unknown  [28]. 
By  analyzing  the  dynamics  transverse  to  the  synchronization  sub-manifold  containing  the  end 
lasers,  sufficient  conditions  for  stability  may  be  derived  in  the  case  of  semiconductor  lasers. 

On  the  other  hand,  in  systems  of  two  delay  coupled  lasers,  lag  synchronization  is  also  ob¬ 
served.  By  averaging  over  the  temporal  dynamics  of  the  outer  lasers  in  our  ring  laser  system,  the 
system  may  be  considered  to  similar  to  a  two-laser  delay  coupled  system,  where  lag  synchro¬ 
nization  is  understood.  That  is,  the  cross  correlation  between  the  two  lasers  peaks  not  at  zero, 
but  at  lag  time  equal  to  the  delay.  This  has  also  been  observed  in  two  laser  systems  [19,21]. 
The  theoretical  analysis  of  this  will  be  presented  elsewhere  for  the  fiber  ring  lasers  in  the  hub 
formation. 
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